NASA-CR-1 96384 


33 “3 


Progress Report 



Submitted to: 
Grant Title: 
Grant Number: 
Organization: 


NASA Langley Research Center 

ROBUST STABILITY OF SECOND-ORDER SYSTEMS 

NAG- 1-1 397 

Georgia Institute of Technology 
Atlanta, GA 30332-0150 


Principal Investigator: Dr. C.-H. Chuang 

School of Aerospace Engineering 
Georgia Institute of Technology 
Atlanta, GA 30332-0150 
(404) 894-3075 

Period Covered: Aug. 24, 1993 to Aug. 24, 1994 

Date of Submission: Aug. 19, 1994 


(NASA-CR-196384) ROBUST STABILITY N94-37590 

OF SECOND-ORDER SYSTEMS Progress 

Report, 24 Aug. 1993 - 24 Aug. 1994 

(Georgia Inst, of Tech.) 23 p Unclas 


G3/63 0019671 



Summary 


This particular research project has a no cost extension until May 23, 1995. The 
extension is due to unexpected funding cut at the second year of this three years (proposed) 
project. This progress report gives current progress of the research in nonlinear robust control 
using positive real concept 

The progress is documented in the following draft paper. In the paper, the manipulator 
dynamics is reformulated differently from the existing equations of motion for free base robots 
This new formulation gives a compact form of the dynamic equations for easy computation. 
The nonlinear terms are now considered. The results show that for an additional nonlinear 
friction term, the feedback controller designed using passivity concept works quite well. 
Although, design of such a controller requires simulation of the dynamics, for the example 
shown in the following draft, this design procedure is feasible. 


NONLINEAR CONTROL OF SPACE MANIPULATORS WITH MODEL UNCERTAINTY 


Manoj Mittal 1 and C. -H. Chuang 2 
School of Aerospace Engineering 
Georgia Institute of Technology 
Atlanta, Georgia 30332-0150 


and 


Jer-Nan Juang 3 
Spacecraft Dynamics Branch 
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Hampton, Virginia 


Abstract 

Nonlinear control using feedback linearization or inverse dynamics for robotic manipulators yields good results 
in the absence of modeling uncertainty. However, modeling uncertainties due to unknown joint friction coefficients 
and payload variations can give rise to undesirable characteristics when these control systems are implemented. In 
this work, it is shown how passivity concepts can be used to supplement the feedback linearization control design 
technique, in order to make it robust with respect to the uncertain effects mentioned above. Results are obtained for 
space manipulators with freely floating base; however, they are applicable to fixed base manipulators as well. The 
controller guarantees asymptotic tracking of the joint variables. Closed-loop simulation results are illustrated for 
planar space manipulators for cases where uncertainty exists in friction modeling and payload inertial parameters. 


1. Introduction 

The dynamics of space manipulators differs from that of fixed base manipulators since their base is free to 
move. The base could be either a spacecraft or a satellite. The movement of manipulator arms produces reaction 
forces and torques on the base. Therefore the resulting motion of the base has to be accounted for in the dynamic 
model for the manipulator. However, Papadopoulos and Dubowsky 1 showed that a dynamic model for space 
manipulators with a free base is similar in structure to the dynamic model for fixed base manipulators. An obvious 
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similarity is that the inertia matrix in each case is symmetric and positive definite. In fact, the dynamic model for 
fixed base manipulators can be viewed as forming a subset of that for space manipulators. 

Some concepts have been proposed for joint trajectory control and inertial end tip motion control of space 
manipulators. Vafa and Dubowsky 2 developed an analytical tool for space manipulators, known as the virtual 
manipulator concept The virtual manipulator is an idealized kinematic chain connecting its base, the virtual base, to 
any point on the real manipulator. This point can be chosen to be the manipulator's end effector, while the virtual 
base is located at the system center of mass, which is fixed in inertial space. As the real manipulator moves, the end 
of the virtual manipulator remains coincident with the selected point on the real manipulator. Additionally, it can be 
shown that the change in orientation in the virtual manipulator's joints is equal to the change in the orientation of the 
real manipulator's joints. While these features give the designer the ability to represent a free floating space 
manipulator by a simpler system whose base is fixed in inertial space, the associated transformation depends on 
knowing the system parameters exactly. Alexander and Cannon 3 showed that the end tip of the space robot can be 
controlled by solving the inverse dynamics that includes motion of the base. Their method assumes the mass of the 
spacecraft to be relatively large compared to that of the manipulator it carries, and also requires much computational 
effort to determine the control input. Note that some future space systems are expected to have the manipulator and 
spacecraft masses of the same order. Umetani and Yoshida 4 proposed the generalized Jacobian matrix that relates 
the end tip velocities to the joint velocities by taking into account the motion of the base. However, robustness of 
the control scheme with respect to modeling uncertainties is not discussed. Masutani et. al. 5 proposed a sensory 
feedback control scheme based on an artificial potential defined in the sensor coordinate frame. This scheme is 
based on proportional feedback of errors in the end tip position and orieniation as well as feedback of joint angular 
velocities. A robust controller based on feedback linearization using a simplified nonlinear model and passivity 
concepts was developed by Chuang, Mittal, and Juang 6 . 

In this paper a robust control scheme using feedback linearization of the complete model and passivity concepts 
is proposed for space manipulators in which joint friction and mass variations are considered. The proposed 
technique can be used for fixed base manipulators also. The control scheme uses inverse dynamics; however, it is 
robust in the face of bounded modeling uncertainties which might arise due to improper or absent friction modeling 
as well as payload variations. The controller asymptotically tracks prescribed time varying joint angle trajectories 

whose acceleration is bounded in the L 2 space. Although only two types of uncertain dynamic effects are 

♦ 

considered hoe, the controller can be applied to systems with other types of uncertainty. 


2. Dynamics of Space Manipulator System with Uncontrolled Base 
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This development of a nonlinear dynamical model for space manipulator systems whose base is uncontrolled is 
discussed in this Section. The development of the expressions for linear and angular momenta of the system closely 
follows that given in Reference [5], however, the form of the final equations of motion is different A space 
manipulator system in the satellite orbit can be approximately considered to be floating in a non-gravitational 
environment As shown in Figure 1, the manipulator and the base can be treated as a set of n+1 rigid bodies 
connected through n joints. The bodies are numbered from zero to n with the base being 0 and the end tip being n. 
Each joint is then numbered accordingly from one to n. The angular displacements of the joints can be represented 
by a joint vector. 


q = [qj q 2 -q„] T 0) 

The mass and inertia tensor of the i^ 1 body are denoted by mj and If, and the inertia tensor is expressed in terms of 
the base frame coordinates. 

2.1 Kinematics 

A coordinate frame fixed to the orbit of the satellite can be considered to be an inertial frame, denoted by Zi. In 
addition to Zi, another coordinate frame Zb is defined that is attached to the base with its origin located at the base 
center of mass. The attitude of the base itself is given by roll, pitch, and yaw angles. In the sequel, all vectors are 
expressed in the base fixed coordinate axes. 

Let Rj and rj be the position vectors of the center of mass of the i* link with respect to frames Zi and Zb, 
respectively. Then 


Ri = Rb + H (2) 

where Rb is the position vector from the base cento - of mass to the origin of the frame Zp Let Vi and Qi be the 
linear and angular velocities of the center of mass of the i* link with respect to frame Zj and vj and a>i be the linear 
and angular velocities of the same point with respect to frame Zb. Then Vj and £2i can be written as 
* 

Vi = Vb + Vi + ne X ri (3) 

fli = Ob + ©i (4) 

Vb and Qg are the linear and angular velocities of the base center of mass with respect to frame Zj. Note that for 
any space manipulator, vj and co; for each link can be represented by the following forms 
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Vi = Ju(q)q 


(5) 


“i =^Ai(q)<i 

where JLj(q) and JAi(q ) e R 3xn are the Jacobian matrices for the i* link. 

The position of the system center of mass with respect to the base frame depends on the joint angles, 
below are two measures related to the system center of mass. 

n 

m c = L m i 
i=0 


r c (q)= 


Invito) 


i=0 




2.2 Linear and Angular Momenta 

The linear momentum P and the angular momentum L of the whole system are defined as follows 

P= tniiVi 

i=0 

L=l[l i £2 i + m i R i xV i ] 

i=0 

Substituting Equations (2) through (8) into Equations (9) and (10) yields 


P = H v V b + H^b + H vq q 


L — hJqV b + H q Q b ■*" Hjjq q + Rb x P 

where * 

H v = m e I, H v e R 3x3 

H yo = - m c[r e x], H vQ eR 3x3 

H vq = im i J u , H vq eR 3xn 

i=l 


( 6 ) 

Given 

(7) 

( 8 ) 

(9) 

( 10 ) 

(ID 

( 12 ) 

(13) 

(14) 

(15) 
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( 16 ) 


T 

Hq = ilj + Im^i x] [rj x], H q e R 3 * 3 

i=0 i=l 

H Q3 = i{l i J Ai +ni i [r i x]J li }, H^eR 3 *" (17) 

i=l 


Also, for any vector 


f = 


ft 

h • 
h. 


(18) 


[f x] is defined as follows 


[ f xj = 



-u 

o 

ft 



(19) 


and I is an identity matrix of appropriate dimensions. 


Since the working environment is non-gravitational and no actuators generating external forces are employed, 
the linear and angular momenta of the whole system are conserved. Since the inertial frame is fixed to the orbit, the 
whole system can be assumed to be stationary with respect to the inertial frame at the initial state. Thus the above 
two momenta are always zero for the whole system. Note that it is implicitly implied that the satellite is a non- 
spinning body. By using the fact then that the linear and angular momenta are zero. Equations (1 1) and (12) result in 


V B =-H; 1 [H va n B + H v<1 q] 

(20) 

= -[Hq - [h*, - Hy^H^jq 

(21) 


2.3 Manipulator Dynamics 

The total kinetic energy of the space robot can be written as 

T = |i(m i V?'v i+ fiTi A ) 

2 i=o 


( 22 ) 


Using Equations (3) through (8) and (13) through (17) the kinetic energy can be expressed as 
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( 23 ) 


where Hq is the inertia matrix corresponding to the fixed base manipulator 

H q = ikjLJu+JLvJ. H q € R nXB (24) 

i=r 

Equation (23) for the system kinetic energy can be simplified as follows. Substituting for Vb from Equation (20) 
leads to 

* T=-in£Mn B + Q^Zq + ±q T Wq (25) 


where 


M = H q -HJ q H; 1 H vQ , MeR 3 ' 3 (26) 

Z = Ho,, - , ZeR 3 " (27) 

W = H q -Hj q H; , H WJ , W e R™ (28) 


Further, substituting for fig fr° m Equation (21), one obtains an expression for the system kinetic energy solely in 
terms of the joint variables 


T = |q T D(q)q, DeR BXB 


(29) 


where D is the inertia matrix of the system and is given by 

D = W-Z t M“ 1 Z (30) 

It can be shown that D = D T > 0. It is interesting to note that the system inertia matrix obtained in Reference [1] is 
of the same form as the above. However, the expressions for W, M, and Z matrices are different This is because a 
different approach, viz. the concept of barycenters, is used in the model derivation of Reference [1]. 
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Since there is no potential energy in non-gravitadonal environment, the Lagrangian A, is equal to the kinetic 


energy 


A = T 


(31) 


So the system dynamics is given by 


dt dq J 9q 


(32) 


where x is an nxl vector of joint torques. The equations of motion for space manipulators can then be written as 


D(q)q + h(q,q) = t 


(33) 


where 


h(q,q) = C(q,q)q + x f (34) 

Paralleling the development for fixed base robots given by Spong and Vidyasagar 7 , the elements of the matrix C are 
obtained as 



dD^ + 

dq { 


d p ij 
Sqj 8q k 




(35) 


and Xf represents the joint torque vector due to friction. As pointed out by Craig 8 , the total friction at each joint can 
be taken to be the sum of Coulomb friction and viscous friction. Coulomb friction is constant except for a sign 
dependence on the joint velocity. Viscous friction, in general, depends on various powers of joint velocity. 
However, higher powers contribute significantly only at high joint velocities. Manipulators usually do not attain 
such high velocities. Therefore, it is sufficient to consider only the linear dependence of viscous friction on joint 
velocity. Figure 2 shows a friction model consisting of Coulomb friction and linear viscous friction. Using this 
model, the joint friction torque vector can be represented as 


x f = 'Fsgn(q) + fq 


(36) 
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where V is a diagonal matrix consisting of Coulomb friction constants for the joints, and T is a diagonal matrix 
consisting of viscous friction coefficients for the manipulator joints. It turns out that in many manipulator joints, 
friction also displays a dependence on joint position. However, such effects are not considered here. There are 
other effects like bending effects that are difficult to model and also neglected in the present model. 

2.4 Base Motion 

The translational velocity of the base center of mass can be written in terms of joint velocities by using the 
expression for Qq from Equation (21) in Equation (20). 

V B =- H; 1 [H vq - H vQ M -1 z]q (37) 


Also, the base angular velocity from Equation (21) is 

n B = -M -1 Zq 


(38) 


Using the above expressions, the evolution of the base position and orientation with time can be determined as 
follows 


V 


CyCg SyC^ CySgC^ SyS^ 

y b 

= 

SyCg SySgS^ "f CyC^ SySgC^, CyS^ 

M. 


-s e CgS* C e c* 


V 


1 s^ tan(0) c^ tan(0) 

8 

= 

0 c^ — s^ 

y 


0 s^ sec(0) c$ sec(0) 


(39) 


(40) 


where c ( .) = cos(-) , s ( ) = sin(-) . 


3. Control System Design 

Assuming that the dynamics of the space manipulator is described by Equation (33), where D and h are 
completely known, the feedback linearization or inverse dynamics 7 technique can be used to design controllers for 
tracking prescribed command trajectories for the joint angles. This can be accomplished as outlined in the following 
sub-section. 
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3.1 Feedback Linearization 


Let the joint torque vector be of the following form. 

't = Du + h (41) 

where u is the pseudo-control, i.e., it is the control input to the resulting linearized system. With the control law 
given by Equation (41), the closed-loop system becomes 



(42) 


where 


A = 


~~o 
1 

D3 

II 

'O' 

LO oj 


I 


(43) 


A simple PD (Proportional-Derivative) type of control law is chosen for the feedback linearized system 


u = q d + K 2 (q d - q) + Kj (q d - q) 


(44) 


where Ki and K 2 are proportional and derivative gain matrices, respectively. These matrices are usually chosen to 
be diagonal in order to achieve decoupled response among the joint angles. Substituting for u from Equation (44) 
into Equation (42), one obtains 


e = A c e (45) 

where e = [e^ e 2 ] T ,ei = q d -q, e 2 =q d -q, A c = A-BK, and K = [Kj K 2 ]. If Ki > 0 and K 2 > 0, the error 

dynamics as given by Equation (45) is asymptotically stable. The freedom in selecting the gain matrices can be 
utilized to meet performance specifications for the closed-loop system. 

The preceding discussion assumes availability of perfect knowledge about the nonlinear system dynamics. 
♦ 

However, in practice, D and h are usually imprecisely known due to modeling inaccuracies. It is assumed that the 
controller is designed using best estimates for friction models and for a nominal value of the end-tip payload. The 
actual joint friction will be different from that assumed in controller design, and the actual tasks might require 
handling a variety of payloads. Thus the controller uses computed versions of D and h. The objective here is to 
design a control law that is robust for bounded variations in D and h due to the uncertain dynamic effects in friction 
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and payload. This issue of robust control design is discussed in the following sub-section, where it will be seen that 
the control law results in closed-loop asymptotic tracking. 


3.2 Robust Feedback Linearization Using Passivity 

In the presence of modeling uncertainties, the control law is given as 

x = D c u + h c (46) 

where D c and hg are computed versions of D and h respectively. Substituting for x and u from Equations (46) and 
(44) into Equation (33) it can be shown that the closed-loop system dynamics is given by 

e = A c e + Bv (47) 


where 


v = Au + 8 


(48) 


and 

A = [l-D~‘D C ], 8 = D -1 [h-h c ] 


(49) 


The first step in the proposed design is to choose the gain matrix K = [Kj 
the linear system given by 


K 2 ] and an output matrix F such that 


e = A c e + Bv 

(50) 

y=Fe 

is SPR (Strictly Positive Real). This can be achieved as outlined in the following theorem. A definition of the 
concept of Strictly Positive Realness can be found in Slotine and Li 9 . 

* 

Theorem 1 [10). Let Ki and K 2 be such that 
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Kj = diag[k H ]; k H > 0, i = 1 n 

K 2 = diag[k 2i ]; k 2i > 0, i = l,...,n (51) 

( k 2i) 2 >k u . i = 


then if F = K, the system described by Equation (50) is SPR. 

Note that the conditions of the Theorem as prescribed by (51) are extremely easy to satisfy. 

With the linear System (50) being SPR, the Passivity Theorem 11 can be used to design asymptotically stable 

(T _ v /2 

controllers as shown in the following theorem. The notation flxlj. = M x T x dt J is used in the sequel. 

Theorem 2. Let the desired trajectory for joint variables be such that q d e L 2 . Further, in the control law given by 
Equation (46), let h c be such that 


D- 1 [h-h c ]L< c|u[ T + d VT>0 (c>0,oo>d>0) 


(52) 


If D c is chosen such that 


where 


D -1 D C 2 arl (a > 0, r > 0) 


(53) 


a>- 


c + 1 


(54) 


then the closed-loop System (47) is asymptotically stable. 

Proof. The closed-loop system as given by Equation (47) can be represented in block diagram form as shown in 
Figure 3. It is first shown that the nonlinear block in the feedback path is passive 1 1 . Consider 

T 

I = J -u T vdt (T > 0) 
o 
T 

= J-u T (Au + 5)dt 
o 
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( 55 ) 


f T „ > 

( ^ 

- f u T Audt 

+ - f u T 6dt 

v 0 > 

V 0 > 


Let the first and second integrals on the right hand side be denoted by Ii and I2 respectively. Then 

Ij = J u T [D~ 1 D e - l]u dt (56) 

0 

Noting that Inequality (53) gives rise to D -1 D C - 1 > (ar - 1)1 , Equation (56) can be used to obtain the following. 


I,*(ar-l)||u|£ 


(57) 


On the other hand. 


-I 2 =|u T D- 1 [h-h c ]dt 

0 

< (julj. |D _1 [h - h c ]L (Holder 1 s Inequality) 

£ |4r [c|ul T + d] (Inequality (52)) (58) 


Hence 


I £ (ar - c - l)|u|£ - d|ul T = f {|u| T ) 


(59) 


It can be shown that if (ar - c - 1) > 0, then 


f (I u 8t ) ^ - 4( _ ■ „ v l-lr ^ 0 m 

which in turn would imply 

T j 2 

j-u T vdt£ - r VT>0 (61) 

0 4(ar-c-l) 


* 

Thus a sufficient condition for the nonlinear block to be passive is that a > (c + l)/r. 

Additionally, the transfer function of the feedforward block [Ac, B, K) is proper and has no poles on the 
imaginary axis. Hence it has finite gain 12 . Since q d e L 2 , then using the Passivity Theorem 11 , one can conclude 
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that the signals u, Ke, and v are bounded. Moreover, since the feedforward block is SPR, Ke(t) = Kje^t) + K^Ct) 
goes to zero asymptotically. This in tum implies that ei(t) and e 2 (t) individually approach zero asymptotically 8 . 

Remark. The conventional method of designing robust controllers for robotics problems involves introducing an 
outer loop control 7 in order to compensate for the uncertain dynamic effects. This has been known to introduce 
chattering in the control. When the design method is modified to make the control smooth, closed-loop asymptotic 
tracking is generally compromised to some extent. In the method proposed in Theorem 2, an outer loop control is 
not employed. Instead, achievement of robustness can be qualitatively understood as follows. Depending on the 
amount of uncertainty contained in h, the computed version of the inertia matrix D, is modified such that conditions 
(53) and (54) of the Theorem are satisfied. 

The results of the Theorem are applicable to space manipulators as well as fixed base manipulators. Finally, it 
should be noted that the control design obtained using the results of the theorem is not unique. Additionally, there is 
a considerable amount of margin for performance optimization. 


4. Simulation Results 

As an example, results of applying Theorems 1 and 2 in order to achieve a robust control design will be 
illustrated for planar space robots. Nonlinear dynamic models for the robots are obtained using the results of Section 
2. The controller is designed for nominal values of joint friction coefficients and payload mass; whereas the actual 
robot will be considered to have variations in these. Simulation is carried out using variable-step fourth and fifth- 
order Runge-Kutta methods. The base and link masses in the following examples are assumed to be of the same 
order of magnitude. Closed-loop results are generated for step command to the joint angles. Note that in general, 
for end tip motion control in the inertial space, the inverse kinematics problem needs to be solved in order to 
generate a command trajectory for the joint angle vector. 

4-1 Planar One Link Space Manipulator 

Figure 4 shows a planar one link space robot. Equation (33) describes the dynamics of this one degree of 
freedom system. The system inertia, computed using Equation (30), turns out to be 

D(q 1 ) = W u -Z| 1 /M 3 . 3 (66) 

where 

W u = m 0lP f + I 1 (67) 
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Z 3 1 - m 01 P! (p 0 Ci + Pi ) + Ii 


( 68 ) 


M 3 3 — m 01 (po + Pi + 2poPiCj ) + 1 0 + 1] (69) 

Using Equations (34) and (35), h is determined to be 

h(Qi ,qi ) = m °'^ 2 Pl ' 1 [ m oiPo(Po + PiCi ) + 1 0 ] ■ [moiPi (Po^i + Pi ) + Ii ]q? + Vsgnfai ) + TCi (70) 

M 33 

where y and y are respectively the coefficients of Coulomb and viscous friction. In Equations (67) through (70) 
m oi = m o m i /m c , m c = mo + mi, Cj = cos(qj), and Sj =sin(q 1 ). The quantities mj, Ii, and pj are defined for an 

“equivalent” link that is obtained by removal of the end dp payload and absorbing its inertial characterisdcs in those 
of the link itself. The concept of the equivalent link is explained in the Appendix. 

It can be seen easily that as m 0 — and I 0 -» 


D->mipf+I lt h->\|/sgn(qi) + yqi (71) 

which represents the case of a fixed base manipulator. Equation (39) is used to determine the evolution of the base 
position with time 


Xi. = 




m. 


P1S V 1 (P() S y + Pl$yl ) 


3,3 


Hi 


Yb=- 


Ei 

m„ 


2 

— Pl C <|>l ■*" ^ (Po^y "*"PlCyl) 


(72) 


Qi 


where s yl = sin(y + q, ), c yl = cos(y + q t ) . Finally, the base attitude dynamics is obtained using Equation (40) 


• ^3,1 • 

v = "m q ‘ 

^3,3 


(73) 


Next, a feedback controller is designed for the one link space manipulator using the results of Theorems 1 and 
2. Dc and h c are obtained from Equations (66) and (70) respectively, by using the nominal values of end-link 
parameters and friction coefficients: 
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D c *p^ 


*>c = H» i-»m k .p, -»Pv..I, -*Ii c .'t'-*V c .7-*Tc 


( 74 ) 


Closed-loop results are generated for a step command of 1 radian to the joint angle. Table I lists physical parameters 
of the example robot used in simulation. The feedback gains are chosen to be kj = 0.4 and k 2 = 1.0. This choice of 
gains satisfies Conditions (51) and in case of no modeling uncertainty, yields a closed-loop response without any 
overshoot Simulation results will be shown for the case in which there is no modeling uncertainty, and for another 
case that involves uncertainty. 

Figures 5 through 8 show closed-loop responses for the nominal case. Figure 5 shows that asymptotic tracking 
in the joint angle is achieved. Figures 7 and 8 show that the base moves in reaction to link motion; this is due to the 
conservation of linear and angular momenta as discussed in Section 2. However, the joint angle still achieves the 
appropriate commanded value. The closed-loop responses for the case involving uncertainty will be overlaid on 
Figures 5 through 8. 


5. Conclusions 

A robust control method based on feedback linearization and passivity concepts is proposed for space 
manipulators. The method is applicable to fixed base manipulators as well. The control law results in asymptotic 
joint angle tracking in the face of bounded uncertainties in friction modeling and payload inertial characteristics. 
For the first time, closed-loop simulation results will be presented using this control method. 
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A ppendix 

The concept of an equivalent link is described here. Following are two equivalent representations of the same 
link. The Figure on the left shows the actual link, while the one on the right shows a hypothetical equivalent link 
that is obtained by removal of the end tip payload and absorbing its inertial characteristics into the link. 
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L 


J 


Actual Link 


Equivalent Link 


The geometric and mass properties of the equivalent link are given as follows: 


nij = nij +m e 


_ m^+m.l, 

Pi “ ' 

nij +m c 


Il = I ;+ -p'j ) 2 

mj +m e ' ’ 
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Linki 


Base Frame 



Earth 


Figure 1 . A Space Robot 
Friction Torque 
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Figure 7. Base Attitude Response for the Single Link Planar Space Robot. 



Figure 8. Motion of the Base Center of Mass for the Single Link Planar Space Robot. 




Table I. Physical Parameters of Single Link Robot 


Body 

p (meter) 

1 (meter) 

m (kg) 

I (kg.m 2 ) 

0 (Base) 

3.0 

- 

5.0 

30.0 

1 (Link) 

3.0 

6.0 

1.0 

3.0 
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